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Abstract 


We  construct  a  family  of  random  step  functions  {x  (t)j  whose  members 

all  have  the  same  power  spectrum  and  such  that  as  ri->oo,  x  (t)  converges 

to  X  (t),  the  Gaussian  process  with  the  same  spectrum,  We  illustrate  the 

procedmre  for  calculating  the  general  multivariate  distribution  of  the 

processes  (x  (t)]  by  calculating  the  univariate,  bivariate  and  trivariate 

distributions.  We  show  how  a  suitably  constructed  univariate  entropy  can 

serve  as  an  index  of  the  extent  to  which  x  (t)  has  approached  the  Gaussian 

n 

limit  x_(t). 

00 
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I,  Introduction 

Consider  a  zero-mean  stationary  random  process  x(t)  with  correlation 
function  CCf ),  i.e.   a  process  x(t)  for  which 

E  x(t)  -  0,    E  x(t)x(t+t')   -  C{T)       , 

where  E  is  the  expectation  operator.     A  considerable  part  of  noise  theory  is 

devoted  to  study  of  the  case  where  x(t)  is  Gaussian.     In  this  case,   the  law 

of  x(t),  i.e.  the  probability  that  x(t   )<c,   l<n<N,  1<N<  oo,  where 

nn—      —  — 

the  c     are  arbitrary  real  numbers,  is  the  familiar  N-varlate  Gaussian  dis- 
tribution whose  moment  matrix  can  be  expressed  simply  in  terras  of  G(t')'-  -^  • 
Thus,   in  the  Gaussian  case,   C(r)  uniquely  specifies  x(t),   which  is,   of  coiirse, 
the  particular  beauty  of  this  case.     More  generally,   G(f)  gives  only  a  more 
or  less  incomplete  characterization  of  x(t). 

The  object  of  the  present  paper  is  to  describe  and  study  an  infinite 
f  amily  •[  X  (t)]  ,   l<n<  oo,   of  non-Gaussian  random  processes,  which  converge 
to  a  Gaussian  process  x     (t)  as  n  ->  oo,  and  whose  laws  are  calculable,   at  least 
in  principle.     The  family  {x  (t)]   will  be  constructed  in  such  a  way  that  all 
its  members  have  the  same  correlation  function  C(f),   or  equivalently,   the 
same  power  spectrum 

00 

J(co)  -  2^     I       exp(-ioof)C(f)dr 

-00 

Briefly,  we  shall  say  that  /x  (t)]  is  an  isospectral  family.  Since  as  n  ->oo, 
the  general  appearance  of  the  sample  functions  of  x  (t)  changes  and  approaches 
that  of  the  Gaussian  process  x  (t)  with  the  same  power  spectrum  (or  correla- 
tion function),  the  general  inadequacy  of  J(co)  or  C(f)  as  a  means  of  character- 
izing random  processes  and  the  general  need  for  higher-order  statistics  are 
put  into  rather  strong  focus. 
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The  method  we  use  to  construct  x  (t)  from  an  underlying  shot  noise 
(Poisson  process)  is  a  familiar  one.  The  novelty  of  ovoc   treatment  consists 
in  observing  that  when  the  shots  are  rectangular  the  members  of  f x  (t)| 
can  be  described  in  theoretically  complete  detail. 


II.  Construction  of  fx  (t)i 


Let  the  sequence  t.  (j"  ...,-1,0,1,...)  be  the  occurrence  times  of  the 
shots  in  a  stationary  shot  noise  (Poisson  process)  with  an  average  rate  of 
p  shots  per  second.  More  specifically,  we  have  the  following  five  properties 
(among  others): 

1.  The  occurrence  of  m  shots  in  the  interval  I  and  the  occurrence  of 
n  shots  in  the  interval  I  (m,n"0,l,2,...)  are  independent  events  if  I  and 
I  do  not  overlap. 

2.  The  probability  of  one  shot  in  an  infinitessimal  interval  of  length 
At  is  pAt  ♦  o(At),  whereas  the  probability  of  more  than  one  shot  is  o(At). 

3.  The  probability  of  m  shots  in  any  interval  of  length  T  is  given  by 

the  Poisson  distribution 

X™   -\ 
p(m/X)  ■  ijj-  e    »  m  -  0,1,2,...  , 

where  X  -  pT,  the  parameter  of  the  Poisson  distribution,  is  the  common  value 

of  the  mean  and  variance  of  a  random  variable  which  takes  the  values  m  ■  0,1,2,..., 

with  probabilities  p(ra;X). 

U.  Given  that  m  shots  have  occurred  in  &n  interval  of  length  T,  their 
occurrence  times  t,,...,t  (without  regard  to  order  of  appearance)  are  inde- 
pendent,  identically  distributed  random  variables  with  the  uniform  probability 
density  l/T. 

5.  The  probability  density  of  the  interval  between  successive  shots  is 
p  exp  ("pt). 
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For  discussion  and  derivation  of  these  properties  we  refer  the  reader 
elsewhere H,  [3],  M, 

Now  let  h(t;a)  bo  a  step  function  of  unit  height  and  width  a,  i.e. 

h(t;a)  -  1  ,  0  <  t  <  o  , 

(1)  "   " 

h(t;a)  ■  0  ,  -co  <t<0,  a<t<oo. 

To  construct  the  isospectral  family  {3e(t)jwe  write 

(2)  x„(t)  -  (1/  Vn)  F      h(t  -  tJ"V)  -V^a,     l<n<oo, 

"^  j-  -co  ^ 

(n) 

where  the  t .       are  random  times  belonging  to  a  shot  noise  with  an  average 

rate  of  n  shots  per  second.     It  is  clear  that  each  x  (t)  is  a  stationary 
random  step  function.     To  calculate  expectations  of  lagged  products  of  x   (t) 
we  follow  the  usual  procedure:     First  we  use  property  U  (above)  to  calculate 
conditional  expectations  for  a  long  finite  interval  of  length  T  containing 
exactly  N  shots  and  then  we  average  over  N  using  property  3.     Again  we  refer 
to  the  literature  for  details'-  -•,     In  particular  we  find 

E  X  (t)  -  0,     E    x^(t)x^(t+r)  -  c(-cO,  1  <  n  <  00  •  , 
where  C(r)  is  given  by 

00 


(3) 


C(f)  -  I     h(tja)h(t+rja)dt  -  a  -   |r|,  0     |r|  <  a  , 


C(r)  -  0  ,   |ti  >  a  . 
The  corresponding  common  power  spectrum  of  the  isospectral  family  j x  (t)  |  is 


00 


a)  -  jj  J  exp(-l(»r)C(r)  -  (2/nca^)  8in^(coa/2)  . 


-00 


-  u- 


The  conatartt  o  is  still  at  onr  disposal.     In  order  to  make  the  members 
of  (x   (t)j  corresponding  to  small  values  of  n  drastically  non-Gaussian,  we  can 
choose  a  «  1,   so  that  there  is  negligible  overlap  of  the  shots  making  up 
X  (t)  if  n  «  l/a»     Then  for  n~l/a  overlap  of  the  shots  begins  to  be 
appreciable  and  as  n  increases  further  the  process  x  (t)  approaches  the 
Gaussian  process  x     (t)  with  the  same  power  spectrum  • 


in.  Typical  Sajiple  Functions  of  (x^(t)} 


To  construct  typical  sample  functions  of  x  (t)  we  first  use  random 
number  tables    -"   to  find  sample  values  of  a  random  variable  4  which  is 
unifomly  distributed  in  the  unit  interval   (0,1).     We  then  observe  that 
the  new  randan  variable 

(k)  n-  ^  log  11^ 

has  the  probability  density  p  exp  (-pt)  of  the  interval  between  successive 
shots  in  a  Poisson  pi\>ces3  with  average  rate  p.     To  see  this,   note  that 

Prob[F"-'-(5)  <  t]    -  Prob[4  <  F(t)]    -  F(t)       , 

so  ttet  f'^(4)  has  the  distribution  function  F(t).     If  F(t)  is  to  be  the 
distribution  function  corresponding  to  the  probability  density  which  is 
p  exp  (-pt)  for  t  >  0  and  zero  otherwise,  then  F(t)  ■  1  -  exp(-pt),  whence 
(U)  follows  at  once.  Moreover,   if  4  is  uniformly  distributed  in  (0,1)  so  is  1-$. 
Thus,  sample  values  of  the  random  variable 

(5)  n„  •  I  log  ? 

generate  time  markers  locating  the  shots  in  a  Poisson  process  with  an  average 
rate  of  p  shots  per  second.     Specifically,   if   "^jji *  •  •  • »  ""?  nm  ^^^  ™  sample 


♦    Specifically,  the  type  of  convergence  we  have  in  mind  is  convergence  in 
distribution,  i.e.  the  law  of  x  (t)  converges  to  that  of  x     (t). 


-  ^- 

values  of  (5),  we  get  m  random  times  tJ^  ,...,t^   by  writing 

Finally,  with  these  values  of  the  random  times  we  use  (2)  to  construct  sample 

functions  of  x  (t). 
n 

Figs.  1,  2  and  3  show  the  results  of  this  procedure  for  the  cases  of 
low,  medium  and  high  density  shot  noise,  respectively.  In  each  case  o  was 
chosen  to  be  0.25*  Fig*  1  shows  a  U  second  sample  of  the  process 

x,(t)  -  r  h(t  -  tj^^;0.25)  -  0.25  , 


Fig.  2  shows  a  U  second  sample  of  the  process 

\^^^   ■  |E   h(t-tj'*^0.25)  -  0.50  , 
and  Fig.  3  shows  a  1  second  sample  (drawn  on  a  different  scale)  of  the  process 

00 


^6^*)  •  J  E   h(t  -  tj^^^-0,25)  -  1.00  . 


The  circle  in  Fig.  2  shows  a  level  where  the  sojourn  of  Xi  (t)  was  so  brief 
that  it  could  not  be  indicated  on  the  scale  of  Fig.  2.  A  similar  reraark 
applies  to  ttie  three  circles  appearing  in  Fig.  3«  In  each  case  the  sample 
functions  were  given  the  initial  value  of  -  ^  o,  i.e.  -a25  for  JE^(t),  -O.^O 
for  Xj  (t)  and  -1.00  for  x,^(t).  When  starting  the  sample  functions  of  X4  (t) 
and  x.At)f   atypical  portions  lasting  about  o  ■  0.25  seconds  occur,  due  to  the 


♦ 


fact  that  there  are  no  shots  before  t  •  0  .   (Theoretically,  each  x  (t)  should 


n 


begin  in  the  infinitely  ronote  past).  To  eliminate  this  transient  behavior,  we 
suppressed  the  first  second  of  the  record  of  x.  (t)  and  the  first  O.U  second  of 
the  record  of  3C-x(t), 


*  This  effect  is  unimportant  for  x, (t),  since  in  this  case  the  overlap  of  shots 
is  sUght. 


x,(t) 
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Figure  1 
Sample  xunction  of  x.  (t)» 
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Figure  2 
Sample  function  of  x,  (t). 
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Figure  3 
Saopl*  fuQotlon  of  x.g(t)« 
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IV.     Probability  Distrlbtttions  of  {Xj^C*)} 

The  family  (x  (t)\  has  the  desirable  feature  that  we  can  calculate  the 
law  of  each  x  (t),   I.e.  all  the  multivariate  probability  distributions  of 
X  (t),  although  the  amount  of  work  required  to  calculate  high-order  probability 
distributions  is  enormous.     We  shall  illustrate  the  general  procedure  by 
deriving  formulas  for  the  univariate,   blvariate  and  trivariate  probability 
distributions  of  x^Ct). 

1.     Univariate  distribution 

By  (1)  and  (2)  the  range  of  possible  values  of  the  random  variable 
X  (t),  t  fixed  ,   is  the  lattice 

T        -  m/  vli  -  v^  a  ,       B  •  0,1,2, ... 
•'m,n  *  »   »   » 

The  probability  that  x  ■  y    is  just  the  probability  that  in  the  a  seconds 

preceding  t  precisely  m  shots  occur.  Consequently,  denoting  the  probability 

that  X  ■  y    by  P  (m),  we  have 
n   •'M,n  ^  n^  * 

/z\  n  /   \         f         \       (not)   -na 

(6)  P„(m)  -  p(m;na)  -  >  /  e 

n  ni4 

where  p(m|na)  is  the  Poisson  distribution  with  parameter  no.  Elementary 

calculations  verify  that 

00  00      A 

as  required.     Note  that  P  (m)  can  be  regarded  as  the  sum  of  n  independent, 
identically  distributed  Poisson  randc«n  variables  with  parameter  a.     Thus, 
the  convergence  of  x  (t)  to  a  zeix)-mean  Gaussian  random  variable  with 
variance  a  is  a  particularly  simple  case  of  the  central  limit  theorem. 
Note  also  that  x  (t)  can  take  arbitrarily  large  positive  values,  but  no 

♦      The  value  of  t  is  irrelevant  since  x  (t)  is  stationary,   so  that  we  can  drop 

the  argument  t. 
♦*     The  fact  that  the  Poisson  distribution  p(jn|X)  approximates  the  Gaussian 

distribution  for  large  values  of  X  is  noted  by  Feller,   0£.  cit.,   p.  176. 
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negative  values  less  than  '  >/n  a.}   this  asymmetry  of  the  sample  functions 
disappears  as  n  ->oo. 

2.  Bivariate  distribution 

The  probability  that  x  (t)  -  y.   while  x  (t+t)  -  y   will  be  denoted  by 
P  (j2,m|T)  and  is  independent  of  t,  since  x  (t)  is  stationary.  The  event 
that  precisely  m  shots  occur  in  the  interval  (u,v)  will  be  denoted  by  E  (u,  v). 
Then  clearly  P  (j^jm^r)  is  the  joint  probability  of  the  events  E  (-o,0)  and 
E  (f  -  a,r).  If  |f|  >  a,  these  two  events  are  independent  and  we  have 

P„(^"*r)  .  P„(£)P„(.)  -  (^  (^  ,-2-  . 

On  the  other  hand,  if  t'-  0,  we  obviously  have 

P^U,m;0)  -  P^a)6,^  , 
where  6^  is  the  lOronecker  delta.  The  interesting  case  is  0  <  |f I  <  o,  for 
then  E  (-a,0)  and  E  (f-  o,r)  are  not  independent  events,  since  the  intervals 
(-a,0)  and  ( tr-  a,T)  overlap.  In  this  case  we  have 

?Jj^,mit)   -  Yi  Prob  |E^,g(-o,r-a),Eg(r-a,0),E^^^(0,^)[ 
if  0  <2'<  a,  and 

PnU»"'i'^)  "  n  Prob  |Ej^^g(t'-a,-a),Eg(-a,t-),E^_g(r,0) 
if  -G  <^<  0|  the  summations  aj-e  over  all  values  of  the  non-negative  integer 
s  compatible  with  the  given  values  of ^ and  m.     In  either  case,  we  find 


"  s-0  (i-s)i  si   (m-s)l 


m-s 


(7) 

.  exp(-n|r|-na)     ,     0  <   |t|  <  a 


-  8  • 
01"  course,  P  (i,ra;'2r)  must  satisfy  the  normalization  condition 

(8)  Z,    P„C^»»;^)  - 1 

for  all  n  andr.     Fo^^'-  0  or  |r|  >  a,   it  is  trivial  to  verify  (8)  directly, 
but  for  0  <   \t\  <  a  it  appears  nrach  harder  to  give  an  independent  proof  of  (8). 


P  (i,  m}f)  must  also  satisfy 

00 


^5^" 


U,B|t-)  -  p^Cb)   . 


I' 

For  small  values  of  m  this  relation  can  be  verified  directly,   but  the  algebra 
becomes  quite  complicated  for  large  ra.     A  more  complicated  identity  st^ns 
from  the  fact  that  C(t)  as  given  by  (3)  must  be  the  correlation  function 
of  every  x  (t),  whence 

Fort'"  0  or   \zr\   >  a  it  is  trivial  to  verify  (9)  directly,  but  for  0  <   \'V\  <  a, 
(9)  becomes  highly  non-transparent.     A  typical  identity  derivable  from  (9)  is 
00       ^^    min(i,m)  gS  ^ 

^^  7^  So      ^-°"  "  ''-°^'  ■  2  '^''^''- 

(Set  n  -  o  -  1,  r-  1/2.) 

3o  Trivariate  distribution 

It  is  clear  by  now  how  the  construction  proceeds  in  general.  Therefore, 

we  calculate  only  the  probability  P  (k,i,m;t','?r« )  of  the  joint  event  that 

x^(t)  -  y.   .  x^(t+t-)  •  To  „  and  x^(t+r+r)  -  y^  .   assuming  that  0  <  r+t'  <  a. 
n      K,  n   n       i: ,  n     n         m,  n 

In  this  case  we  have 

P^(k,j?,m|r,r')  -  Prob{Ej^(-a,0),E^(r-a,r),Ejjj(r-Klr'-o,r+t-')} 
-  2IP^ob[EJ^_J._g(-a,r-a),E^(tr-a,r+t•-a),E^(1^+t^•-a,0), 
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where  we  have  gone  over  to  non-overlapping  events,   and  the  summation  is 
over  all  values  of  the  non-negative  integers  r  and  s  compatible  with  the 
given  values  of  k,  Z  and  m.     Specifically,  we  find 

Bin(k:.s,^-s)        min(k,£,«)  (^)lc-r-s  („^,)r  (^(^.^.^))8 

.  .i — L .i i exp(-nr-nt»-na)       • 

(je-r-s)i       (m-i*r)l 

One  can  use  the  fact  that  P   (k,i,m;r,  t»)  is  a  probability  distribution  to 
generate  a  series  of  very-  complicated  identities,  e.g. 

X:        P  (k.^.ra^t.r^)  -  1, 
k,j2,m-0 

£    P  (k,^,m;2',f)  -  P  (k,/;^)     , 
etc. 

The  9xa»ples  given  make  it  clear  how  to  calculate  P  (k,J?,ra;'t,'er')  for 
other  ranges  of  the  parameters  '^,  t'  and  how  to  calculate  the  general 
N-variate  distribution  P  (k,>««»»^gf  ^f  ••»t|j^^)«     Of  course,   for  large  Nj 
calculation  of  the  N-variate  distribution  is  very  formidable  and  would  require 
the  services  of  a  high  speed  electronic  computer. 


V.    Entropy  of  [x^(t)j 


The  fact  that  as  n  ->oo,  x  (t)  converges  in  distribution  to  the  Gaussian 
process  x     (t)  with  the  same  power  spectrum  follows  from  the  work  of  Fortet'-  -• 

An  index  of  the  closeness  of  x  (t)  to  x„(t)  is  furnished  by  the  entropy  of 

n       00 

X  (t)|  for  simplicity  we  consider  only  the  univariate  entropy  of  x  (t). 
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Since  the  entropy  of  the  Gaussian  random  variable  x      with  mean  zero  and 
variance  o  as  defined  in  the  usual  way  is 


"cot") 


-00 


we  must  insist  that  the  univariate  entropy  of  x    converge  to  H     (a)  as 
n->oo.         As  we  have  seen,   the  univariate  distribution  of  x     is  given  by 

P  (m)  5  Prob  f  x  -y      -  n/ v^  -  >/n  ov  -  p(mina)  s  ^"^  ■    e""""  • 

n  \.    ^      KfU  /  m* 

Thus,  at  first  it  rai^t  seem  that  the  appropriate  entropy  to  consider  is  the 
entropy 

00 

(10)  ^n^**^  ■  ~  n    p(in;no)logp(m;na) 

m-0 

of  the  discrete  random  variable  x  •     However,  H^(a)  as  defined  by  (10)  is 

n  '     n 

independent  of  a,   since  a  is  involved  only  in  the  combination  na,  and 

therefore  H  (a)  cannot  converge  to  H_(a)  -  log  y/Sftea.     In  fact,   H  (a)  is 
n  oo  n 

actually  logarithmically  divergent,    as  the  following  simple  qualitative 
argument  shows:     As  already  noted,   as  n->oo  the  randcHn  variable  x     converges 
in  distribvttion  to  the  Gaussian  random  variable  x_^  with  mean  zero  and 

00 

variance  a.     It  follows  that  for  large  n  most  of  the  distribution  of  x 

n 

is  concentrated  in  the  interval  (  -  y/o,  +  ^a),   i.e.  within  one  standard 
deviation  of  zero,  and  that  the  values  of  p(mjaa)  are  approximately  equal  in 
this  interval.  Since  the  possible  values  of  x  are  a  distance  l/v^  apart, 

*  C.E.  Shannon,  "  The  mathematical  theory  of  communication"  ,  jreprinted  in 
the  book  of  the  same  title  by  C.E.  Shannon  and  W.  Weaver,  Univ.  of  Illinois 
Press,  Urbana,  p.  5U;  19U9.  Shannon's  discussion  of  the  difference  between 
the  entropy  of  a  discrete  random  variable  and  that  of  a  continuous  random 
variable  is  highly  relevant  to  the  present  analysis, 

**  Since  3C  (t)  and  x  (t)  are  stationary  we  need  not  retain  the  argument  t  in 

discussing  univariate  quantities. 
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there  are  approximately  2  •/na   values  of  p(mjna),  all  approximately  equal 
to  l/2yna.  The  resulting  estimate  of  the  entropy  H  (a)  of  the  discrete  dis- 
tribution p(m;na)  is 

__   1        1         __ 
-  2  /na  log  "  log  2  /na  , 

2/na     2v(^ 
which  diverges  logarithmically  as  n->oo» 

The  appearance  of  this  difficulty  is  not  surprising  since  the  entropy  of 

00 

a  continuous  random  variable  is  defined  as  -  J  p(x)logp(x)dx,  where  p(x)  is 

-00 

the  probability  density  of  the  random  variable,  whereas  the  discrete  random 

variable  x  used  to  define  H  (a)  has  no  probability  density,  in  spite  of  the 

fact  that  X  converges  in  distribution  to  the  continuous  random  variable  x  • 
n  00 

What  has  to  be  done  is  to  replace  the  random  variable  x  by  a  related 

continuous  random  variable  x  which  converges  to  x   in  density  as  well  as 

(*)  ♦ 

in  distribution  .  To  achieve  this  we  define  x  as  the  continuous  random 

n 

variable  with  the  probability  density 

00 

(11)  p„(x)  -  r"    v/H  p(Mfna)  h(x  -  y^     xl/^S)       , 

where  h  is  the  step  function  defined  by  (1)  and  as  usual  y    ■  m/ \/n  -  \/na» 

In  other  words,  we  replace  the  discrete  random  variable  x  by  the  continuous 

It 

random  variable  x    where  the   "  mass"  formerly  concentrated  at  the  point  y 

n  J  t-  •'in,n 

is  now  uniformly  distributed  ovei^  the  interval   (y       ,y    _.,  );  of  course,   this 
requires  adjustment  of  the  step  height,   which  accounts  for  the  factor  \/n  in  (11). 

{*)  As  usual,  we  say  that  the  random  variable  x     converges  in  (distribution  as 

n->oo  to  the  random  variable  x       if  F  ,  the  distribution  function  of  x^ 

oo  n'  n 

converges  as  n->oo  to  F     ,  the  distribution  function  of  x     ,  at  every 
continuity  point  of  the  latter.     If  x     and  x       have  probability  densities, 
we  say  that  x^  converges  in  density  to  x      if  p  ,  the  density  of  x 
converges  to  p     ,  the  density  of  x     ,  at  every  continuity  point  of  the 
latter.     Of  course,  when  the  density  p(x)  exists,   the  distribution  function 
F(x)  is  given  by 


F(x)  -      i   p(z)dz 


X 

00 
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* 
Clearly,  it  follows  from  the  central  limit  theorem  that  as  n->oo,  x  converges 

in  density  to  the  Gaussian  random  variable  x  .  The  entropy  of  x  is  defined 

in  the  usual  way  as 

00 


H*(a)  -  -  j  Pj^(x)log  p^(x)  dx 


-00 

Substituting  for  p  (x)  from  (11)  we  find  that 

H  (a)  -  -  23  p(in;na)  log  p(mjna)  -  log  >/n  , 
^  m«0 

80  that  the  logarithmic  divergence  of  H  (a)  is  cancelled  out  and  the  depen- 
dence on  a  is  restored.  Since  the  entropy  of  a  continuous  random  variable 
is  a  continuous  functional  of  the  probability  density,  we  must  have 


H  (a)  ->  H_(a)  -  log  /2nea 
n       00 


as  n  ->oo. 


In  Table  1  we  show  the  results  of  calculating  H  (a)  for  the  value 
a  -  0.25  used  in  constructing  Figs.  1,2  and  3.   (All  logarithms  are  to  the 

base  e.)  The  table  shows  quite  clearly  that  as  n  increases  the  entropy  of 

♦  (♦) 

X     rapidly  approaches  the  entropy  of  the  limiting  Gaussian  random  variable  x     • 

TABLE  r 

n  H*(0.25) 

1  0.617 

2  .581 
U  .612 
8  .66U 

16  .700 

32  .71U 

hOO^**^^  .722 

.726 


CO 


(♦)  Note,   hbwever,  the  initial  decrease  of  H^(Oo25). 

(•«*)The  numerical  work  was  done  by  Miss  P.A.  Smith,  using  E.G.  Molina's  tables, 
"  Poisson's  Exponential  Binomial  Limit"  ,   D.Van  Nostrand  Co.,   Inc.;  19U2. 

(■wt*)'rhe  value  1*00  corresponds  to  the  largest  value  of  p(m;na)  available  in 
Molina's  tables. 
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The  univariate  entropies  of  the  processes  x.(t),   Xi  (t)  and  x-/.(t) 
represented  in  Figs.   1,2  and  3  are  all  close  to  the  limiting  value  of  0,726, 
despite  the  fact  that  these  cases  correspond  to  low,   medium  ar*!  high  density 
shot  noise,  respectively.     However,   the  univariate  entropy  still  seems  to  be 
a  useful  supplement  to  the  power  spectrum  as  an  index  of  the  structure  of  non- 
Gaussian  random  processes.      (It  will  be  recalled  that  all  information  about  the 
relative  phases  of  the  harnoraic  components  of  a  process  are  suppressed  in  its 
power  spectrum,  which  limits  the  utility  of  the  power  spectrum  as  a  means  of 
characterizing  non-Gaussian  processes.) 

VI.     Conclusion 

Every  noise  theorist  is  familiar  with  the  way  in  which  dense  shot  noise 
approximates  a  Gaussian  process.     We  have  given  concrete  form  to  this  model 
by  introducing  an  isospectral  family  |x(tH  of  non-Gaussian  processes 
corresponding  to  progressively  denser  shot  noise.     The  merit  of  our  family 
/x  (t)*!  is  that  we  know  how  to  calculate  in  detail  the  statistical  properties 
of  its  raanbers.     We  have  shown  .how  to  do  this  in  cases  where  the  amount  of 
work  required  is  not  prohibitive. 
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Chicago  2b,  Illinois 

Attn:  La Verne  LaGioia,  Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Microwave  Laboratory 
Building  12,  Room  2617 
Culver  City,  California 
Attn;  M.  D.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
75  Cambridge  Parkway 
Cambridge,  Mass. 
Attn:  Mrs.  Lois  Seulowltz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  33,  New  York 

Attn:  Mr.  C.F.  HcElwaln,  Crfjneral  Manager 

International  Business  Machines  Corp. 
Military  Products  Division 
Owego,  Mew  York 

Attn:  Mr.  D.  I.  Marr,  Librarian 
Department  b59 

International  Resistance  Company 
UOl  N.  Broad  Street 
Philadelphia  8,  Pa. 
Attn:  Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,  N.  W, 
Washington  7,  D.  C. 
Attn!  Mr.  Delmer  C.  Ports 

Dr.  Henry  Jasik,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Electromagnetic  Research  Corporation 
711  Ibth  St feet,  N.  W. 
Washington  5,  D.  C. 
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Lockheed  Aircraft  Corporation 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-75,    Plant  A-1,    Bldg.   63-1 

Bnrbank,   California 

Attni     N.  C.  Harnois 

The  Martin  Company 

P.   0.  Box  179 

Denver  1,  Colorado 

Attn:  Mr.  Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attni  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attn:  Mr.  H.  Warren  Cooper 

Mathenatlcal  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 

The  W.  L.  Maxecn  Corporation 

h60  West  3Uth  Street 

New  York,  N.  T. 

Attn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,  St.  Louis  3,  Missouri 

Attn:  R.  D.  Detrich,  Engineering  Library 

McMillan  Laboratory,  Inc. 

Brownville  Avenue 

Ipswich,  Massachusetts 

Attn:  Security  Officer,  Document  Room 


Melpar,  Inc. 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

Attn:  Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,  Massachusetts 
Attn:  N.  Tucker,  General  Manager 
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Pickard  and  Bums,    Inc. 

2U0  Highland  Avenue 

Needhan  9U,   Mass. 

Attn:     Dr.    J.  T.   DeBettencourt 

Polytechnic  Research  and  Development 

Company,    Inc. 

202  Tillary  Street 

Brooklyn  1,   New  York 

Attn:     Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,    Mass. 
Attn:     Dr.   John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:  P.  S,  Carter,  Lab.  Library 

RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attni  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

AdiB-inlatrator 

The  Ramo-Wooldridge  Corporation 
P.O.  Box  b?U?3  Airport  Station 
Los  Angeles  Ii5,    California 
Attn:     Margaret  C.  Whitnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,   Marvland 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Librarv  -  unclassified 

Sage  Laboratories 
1$9  Linden  Street 
Wellealey  81,  Mass, 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Ubrary 

Sandia  Corporation,  Sandia  Base 

P.O.  Box  5800,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Canpany 

Great  Meek,  Long  Island,  New  York 

Attn:     Florence  W.  Tumbull,   Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park.   California 

Attni     Library,  Engineering  Division 

Sylvania  Electric  Products,    Inc. 
7.00  thirst  Avenue 
Waltham  5li,    Mass. 

/Ittn:     Charles  A.  Thornhill,    Report  Librarian 
Walthan  Laboratories  Library 

Systems  Laboratories  Corporation 
lL852   Ventura  Boulevard 
Sherman  Oaks,    California 
Attn:     Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A,  S.  Thomas,  Inc. 

l6l  Devonshire  Street 

Boston  10,  Mass, 

Attn:     A.   S.   Thomas,    President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 


Microwave  Radiation  Company  Inc. 

19223  South  Hamilton  Street 

Gardena,   California 

Attn:     Mr.   Morris  J.  Ehrlich,    President 

Chance  Vought  Aircraft,    Inc. 
93lh  West  Jefferson  Street 
Dallas,   Texas 
Actrii     Mr.   H.   S.  White,    Librarian 

Northrop  Aircraft,    Inc. 
Hawthorne,    California 

Attn:     Mr,   E.   A.   Freltas,   Ubrary  Dept31b5 
1001  E.  Broadway 

Remington  Rand  Univ.   -  Division  of  Sperrv 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,    Pennsylvania 
Attn:     Mr.    John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,    Inc. 
1221h  Lakewood  Boulevard 
Downey,   California 
Attn:     Engineering  Library  U95-115 

North  American  Aviation,    Inc. 
Los  Angeles  International  Airport 
Los  Angeles  Ii5,   California 
Attn:     Engineering  Technical  File 

Page  Communications  Engineers,    Inc. 
710  Fourteenth  Street,   Northwest 
Washington  5,    D.   C. 
Attn:     Librarian 

Philco  Corporation  Research  Division 

Branch  Library 

Ii700  Wlssachickon  Avenue 

Philadelphia  UU,    Pa. 

Attn:  Mrs.  Dorothy  S.  Collins 


Director,  DSAF  Project  RAND 

Via:  Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

Rantec  Corporation 

Calabasas,  California 

Attn:   Grace  Keener,  Office  Manager 

Raytheon  Manufacturing  Company 
Missile  Systems  Division 
Bedford,  Mass. 
Attn:  Mr.  Irving  Goldstein 

Ravtheon  Manufacturing  Company 
Wayland  Laboratory,  State  Road 
Wayland,  Mass, 
Attn:  Mr.  Robert  Borts 

Raytheon  Manufacturing  Comcany 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  G.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  N.  Y, 
Attn:  Engineering  Library 

Thru:  Air  Force  Plant  Representative 
Republic  Aviation  Corp. 
Farmingdale,  Long  Island,  N.Y. 

Rheem  Manufacturing  Company 
9236  East  Hall  Road 
Downey,  California 
Attn:  J.  C.  Joerger 

Trans-Tech,  Inc. 
P.  0.  Box  3U6 
Frederick,  Maryland 


Chu  Associates 
P.  0.  Box  3B7 
Whitcomb  Avenue 
Littleton,  Mass. 

^!lcrowave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Companv 
Missilp  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50b,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  COTmand 

llith  Street  and  Constitution  Avenue 

Washington,   D.   C. 

Attn:     Mr.  Ottlng,   SRT  ^'.. 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Tnt'l  Airport  Box  7U6 
Baltimore  3,   Maryland 
Attn:      Engineering  Library 
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Wheeler  Laboratoriea,  Inc. 
122  Cutter  Mill  Road 
Great  Neck,  New  York 
Attn:  Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Co. 
Box  91 

fiardena,  Callforni* 
Attn;  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  University  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  h,   Callfomia 

Attn:  Dr.  Samuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.  Electronics 
Research  Lab. 

University  of  California 
Electronics  Research  Lab. 
332  Cory  Hall 
Berkeley  h,   California 
Attn:  J,  R,  Whinnery 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
hROO  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr.  I.  E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,  Pennsylvania 
Attn:  Prof.  A.  E.  Heins 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  0.  C.  Dalman 

Tniverslty  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Prof.  M.  H.  Latour,  Library 

Library 

Georgia  Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,    Georgia 

Attn:     Mrs.    J.H.   Crosland,    Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,   303A  Pierce  Hall 

Oxford  Street,   Cambridge  3",   Mass. 

Attn:     Mrs.   E.L.  Hufschmidt,   Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,    Mass. 
Attn:     Dr.  Fred  L.  Whipple 

University  of   Illinois 
Documents  Division  Library 
Urbana,    Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,    Illinois 

Attn:     Dr.   P.  E.   Moyes,   Department  of 
Electrical  Engineering 

T^.  Johns  Hopkins  University 
Homewood  Campus 
Department  of   Physics 
Baltimore  18,  Maryland 
Attn:     Dr.  Donald  E.  Kerr 

Sandla  Corporation 
Attn:     Organization  lli23 
Sandia  Base 
AlbuTueroue,   New  Mexico 


Applied  Physics  Laboratory 
The  Johfls Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,   Maryland 
Attn:     Mr.   George  L.  Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 

Cambridge  39,  Massachusetts 

Attn:   John  H.  Hewitt 

Hassachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J,  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

University  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilantl,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  lU,  Minnesota 

Attn:  Mr.  Robert  H.  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E.  Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tice 

Dept.  of  Elec.  Engineering 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attn:  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

?5  Johnson  Street 

Brooklyn,  New  York 

Attn;  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A.  E.  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhai?dt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:  Prof.  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Mingins 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn;  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn;  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,    Rhode  Island 
Attn:     Dr.  C.  M.   Angulo 

Case   Institirte  of  Technology 

Cleveland,   Ohio 

Attn:     Prof.   S.   Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momlngside  Heights 

New  York,   N.  Y. 

Attn;  Dr.  Schlesinger 

McGlll  University 
4Iontreal,  Canada 
Attn:  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schultz 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsvlvanla 

Institute  of  Cooperative  Research 

3(iOO  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn;  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:   Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  li,  Pennsylvania 
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Ionosphere  Research  Laooratory 

Pennsvlvanla  State  College 

State  College,    Pennsylvania 

/iTTNj     Professor  A.   H.  Waynlck,   Director 


Institute  of  Hathenatlcal  Sciences 
25  Waverly  Place 
New  York  3,   "ev  York 
ATTN  I     Librarian 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monies,   California 
ATTHt     Dr.   Robert  Kalaba 

National  Bureau  of  Standards 

Washington,    D.  C. 

ATTN:     Dr.  W.   K.   Saunders 

Applied  Mathematics  and  Statistirs  Lab. 

Stanford  University 

Stanford,     dllfomia 

ATTN:     Dr.   Albert  H.  Bowker 

Department  of  Physics  and  Astronony 
Michigan  State  College 
East  Lansing,   Michigan 
ATTNi      Dr.    A.    Leltner 

Dni»erslt"  of  Tennessee 
Knoxvllle,   Tennessee 
ATTN:     Dr.   Fre.1  A.  Flcken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN:     Professor  B.H.  Blsslnger 

General  Atomic 

P.  0.   Box  608 

San  Diego  12     California 

Attn:      Dr.    Edward  Oerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,   Mass. 
ATTNi      Dr.   Arnold   Arons 

California   Institute  of     Technology 
1201  E.  California  Street 
Pasadena,   California 
ATTN:     Dr.   A.  Erdelyl 

Mathematics  Deportment 

Stanford  University 
Stanford,  Califnmia 
.A'TN:     Dr.   Harold  Levine 

Unlversltv  of  Minnesota 

Mlnnpapolis  I'l,    Minnesota 

ATTN:      Professor  Paul  C.   Rosenblooei 

Departnent  of  Mathenatics 

Stanford  University 

Stanford,   California 

ATTN:     Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  .Tohr^  Hopkins  University 
8621  Creorgia  Avenue 
Silver  Spring.    Maryland 
ATTN:     Dr.  B.  S.  Gourary 

{?)Exchange  and  Gift  Division 
The  Library  of  Congress 
Vashlngton  25,    D.  C. 

Electrical  Engineering  Pepartnent 
Massachusetts  Institute  of  TecHnolo^ 
Cambridge  39,    Mass. 
ATTN:     Dr.  L.   J.'  Chu 

Nuclear  Development  Associates,    Inc. 
5  New  Street 

'/(hite  Plains,  New  York 
ATTN:      Library 

California   Institute  of  Technology 
Electrical  Engineering 
Pasadena,   Callfcnla 
ATTN:     Dr.   Zohrab  A.   Kaprlelian 
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Dr.   Rodman  Doll 
311  W.  Cross  Street 
Ypsilanti,   Michigan 

California  Inst,   of  Technology 
Pasadena,   California 
ATTN:     Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographlc  Institute 

Woods  Hole,   Mass. 

National  Bureau  of  Standards 
Boulder,   Colorado 
ATTN:     Dr.   R.  Oallet 

Dr.  Solomon  L.  Schtfebel 

3689  Louis  Road 

Palo  Alto,     ©alifomia 

Oniversity  of  Minnesota 
The  Oniversity  of  Library 
Minneapolis  li,   Minnesota 
ATTN:     Exchange  Division 

Department  of  Mathematics 
Oniversity  of  California 
Berkeley,   California 
ATTN:     Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Box  73 

Lexington  73t   Massachusetts 

ATTN:     Dr.  Shon  Chin  Wang,   Roan  C-351 

Helpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kellehar,  Section  Head 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  0.  Hansoon  Field 

Bedford,  Mass. 

ATTN:  Mr.  Francis  J.  Zncker,  CSBD 

Hq.  Air  Force  Caabrldge  Researah  Center 

Laurence  0.  Hanacoa  FiaU 

Bedford,  Haaa. 

ATTN:  Dr.  Philip  Newman,  CRRK 

Mr.  N.  C.  Oerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Developnent  Corp. 
2liOO  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  23? 

II45  Palisade  Street,  Dobbs  Ferry,  N.  T. 
ATTN:  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
Oniversity  of  Maryland 
College  Park,  Maryland 
ATTN:   Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5.  Mo. 

ATTN:   Professor  J.  Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,  California 

Rrandeia  !Iniversity 
Waltham,  Mass. 
ATTN:  Library 


General  Electric  Company 
Mlctowave  Laboratory 
Electronics  Division 
Stanford  Industrial  Park 
Palo  Alto,  California 
ATTN:  Library 


Saiyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,  California 
ATTN:  Dr.  John  B.  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  California 
>TTN:  DrfladrgoB  Q,  Weill 

Naval  Research  Laboratory 

Washington  25,  D.  C. 

ATTN:  Henry  J.  Passerlni,  Code  5278A 

Dr.  George  Rear 
5  Culver  Court 
Orlnda,  California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,  New  York 

ATTN:  Dr.  Nathan  Marcuvlti 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
85  Livingston  Street 
Brooklyn,  New  York 
ATTN;  Dr.  Jerry  Shmoys 

Department  of  Mathematics 
Oniversity  of  New  Mexico 
Albuquerque,  New  Mexico 
ATTN:  Dr.  I.  Kolodnor 

Mathenatics  Department 
Polytechnic  Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,  New  York 
ATTNi  Dr.  Harry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Grounds 
Aberdeen,  Maryland 
ATTN:  Dr.  Pullon  Keats 

Dr.  Lester  Kraus 
b935  Whitehaven  Way 
San  Diego,  California 

Oniversity  of  Minnesota 
Institute  of  Technology 
Hlnneapclis,  Minnesota 
Attn:  Dean  Athelston  Spllhaus 

Ohio  State  University 
Columbus,  Ohio 
Attn:  Prof.  C.  T.  Tai 
Department  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,  D.  C. 

Attn:  W.  S.  Aaent,  Code  5271 

Naval  Research  laboratory 
Washington  25.  D.  C« 
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